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Introduction: Reading is an acquired-developmental ability that relies on intact language 
and executive function skills. Verbal fluency tasks (such as verb generation) also engage 
language and executive function skills. Performance of such tasks matures with normal 
language development, and is independent of reading proficiency. In this longitudinal fMRI 
study, we aim to examine the association between maturation of neural-circuits support- 
ing both executive functions and language (assessed using verb generation) with reading 
proficiency achieved in adolescence with a focus on left-lateralization typical for language 
proficiency. 

Methods: Normalized fMRI data from the verb generation task was collected from 16 
healthy children at ages 7 11, and 17years and was correlated with reading scores 
at 17 years of age. Lateralization indices were calculated in key language, reading, and 
executive function-related regions in all age groups. 

Results: Typical development was associated with (i) increasingly left-lateralized patterns 
in language regions (ii) more profound left-lateralized activation for reading and executive 
function-related regions when correlating with reading scores, (iii) greater involvement of 
frontal and parietal regions (in older children), and of the anterior frontal cortex (in younger 
children). 

Conclusion: We suggest that reading and verb generation share mutual neural-circuits dur- 
ing development with major reliance on regions related to executive functions and reading. 
The results are discussed in the context of the dual-networks architecture model. 
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INTRODUCTION 

The normal course of development involves continued improve- 
ment in higher-order cognitive abilities (e.g., executive func- 
tions), that originate in prefrontal cortex (PFC) (Diamond, 2002 
for review). There is strong evidence supporting the continued 
anatomical and functional maturation of the PFC related to devel- 
opment of these skills (Sowell et al., 1999). The cognitive skills 
gathered under the umbrella term "Executive functions" include a 
wide range of conscious processes involved in monitoring and 
optimizing performance (Baddeley, 1986) such as error mon- 
itoring, selective attention to the relevant stimulus, inhibiting 
execution of undesired responses, working memory (WM) abili- 
ties (the ability to manipulate and maintain information) , shifting, 
and fluency. 

Fluency is part of the information processing component of 
executive functions and reflects the integrity of neural connections 
and the functional integration of frontal systems. Fluency can be 
evaluated by the speed, quantity, and quality of output (Ander- 
son, 2002). Semantic verbal fluency ability is typically assessed 
either by semantic category fluency (generating members of a cat- 
egory) or by verb generation in response to auditory or visual 



nouns. Individuals performing a semantic category fluency task, 
exhibit activation of the left hemisphere regions related to lan- 
guage and articulation [inferior frontal gyrus (IFG), representing 
word retrieval; middle frontal gyrus (MFC) attributed to verbal 
WM, and supplementary motor area (SMA), which reflects atten- 
tion and motor planning] (Paulesu et al, 1997; Hugdahl et al, 
1999; Pihlajamaki et al., 2000; Gaillard et al., 2003) in both chil- 
dren and adults (Gaillard et al., 2003). The verb generation task 
also reflects semantic verbal fluency ability (Piatt et al., 1999). In 
this task, participants are presented with a noun, either audito- 
rily or visually, and are asked to generate related verbs. Similar 
to semantic category fluency, the verb generation task has been 
found to activate several frontal regions [left and right IFG (BA 
44, 45, 46, 47), and the left medial frontal gyrus (MFG) (BA 6, 
8, 9)] as well as language regions in the brain [left and right 
medial temporal gyrus (BA 19/39, 21), right and left inferior 
temporal gyrus (ITG) (BA 19/37), and superior temporal gyrus 
(STG) (BA 22)] in participants aged 5-18 years of age (Hol- 
land et al, 2001, 2007; Karunanayaka et al, 2010, 2011). We 
have also documented developmental changes in lateralization 
of the language-network supporting the verb generation task in 
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both a cross-sectional (Holland et al., 2007) and a longitudinal 
study (Szaflarski et al., 2006). Specifically, the network becomes 
increasingly left-lateralized with age. Karunanayaka examined 
developmental changes in neurocognitive modules underlying 
the verb generation task using Independent Component Analysis 
(Karunanayaka et al., 2011). Two primary neurocognitive mod- 
ules were identified in that work: (1) word processing [audi- 
tory (BA 22), and then phonological processing (BA 45/47, BA 
44/45/46/9, and 39/40)] and; (2) word generation [the phonologi- 
cal information becomes associated with a semantic meaning (BA 
30/35/19/27/39) a process that can engage visual imagery (BA 17)] 
(Karunanayaka et al., 2011). Although these regions were found 
to be active in the verb generation task from 5 to 18 years old, 
independent component analysis, and structural equation model- 
ing revealed that the connections between STG and the IFG and 
within the frontal lobe itself (IFG and MFG) increase and become 
more left-lateralized from 5 to 18 years. This is in line with the 
timecourse of continued PFC maturation (Holland et al, 2001; 
Karunanayaka et al., 2011), and consistent with studies showing 
increased left-lateralization for verb generation from the age of 5 
to 18 years (Holland et al., 2001). Utilizing this same large imaging 
and behavioral dataset, we can build on the previously published 
results from this task to explore the relationship between this 
developing circuitry underlying verb generation and reading skills. 

Reading is also a linguistic ability, though unlike the phono- 
logical and semantic skills that underlie verbal fluency tasks, it 
must be explicitly acquired. Since reading is a phylogenetically 
new skill (roughly, 5000 years old), neuroscientists suggest that 
the brain uses regions, which are originally devoted to other cog- 
nitive functions and utilize them for gaining information from 
written symbols (see also Price, 2012; Vogel et al., 2013, 2014). 
Reading demands the translation of abstract graphemes (letters) 
into their corresponding spoken language sounds (phonemes) 
that format meaningful words (semantics). It therefore engages 
visual (cuneus and fusiform gyrus) (McCandliss et al., 2003), audi- 
tory/phonological (angular gyrus, STG) (Dehaene-Lambertz et al., 
2002), and semantic brain regions (IFG) (Newman and Joanisse, 
2011) mainly in the left hemisphere (Chiarello et al, 2009). An 
efficient synchronization of these processes is needed for fluent 
reading, which demands the engagement of the executive system 
(Breznitz, 2006). Behaviorally, reading performance has been sug- 
gested to correlate with executive function skills (Christopher et al., 
2012; Kieffer et al, 2013; Booth et al, 2014). Many studies have 
focused on specific aspects of executive functions expected to be 
related to reading. For example, attention shifting and inhibition 
(Kieffer et al, 2013), WM, and speed of processing (Christopher 
et al, 2012), as well as verbal fluency (or verb generation ability), 
all were found to be correlated with reading ability (Snowling et al., 
1997). These findings suggest a major role for executive functions 
in the reading process. 

The neural architecture supporting executive function is 
described by the dual-networks top-down model (Dosenbach 
et al., 2008). This model proposed two cognitive control/ executive 
functions networks with different neuroanatomical correlates. The 
first is the rapid adaptive control network, which is in charge 
of allocating attention to a cue; uses feedback to affect process- 
ing and involves a frontal-parietal circuit. The second network 



is the set-maintenance network that maintains task goals, which 
involves a cingulo-opercular circuit. The activation of these net- 
works changes with development due to changes in connectivity 
between key elements within the network (Dosenbach et al., 2008) . 
Both networks are engaged during reading (Ihnen et al., 2013), but 
only the functional connectivity between the fronto-parietal net- 
work and the fusiform gyrus was found to be positively correlated 
with reading skill and age (see Vogel et al., 2014). 

In this current study, we aim to examine the relations between 
the maturation of neural networks supporting executive func- 
tions and language (during the verb generation task) and reading 
skill. We focus on the language-network described above (Hol- 
land et al, 2007; Karunanayaka et al, 2010, 2011), that includes 
reading regions related to word recognition (i.e., word form area, 
BA 37), and regions related to semantic WM (dorsolateral PFC; 
BA 9, anterior PFC; BA 10), which are part of the executive func- 
tion network in the dual-network model (Dosenbach et al., 2008). 
To explore the role of early emerging language networks (Ahmad 
et al., 2003) in support of later acquired reading skills, we employed 
a novel longitudinal data analysis model that examines the associ- 
ation between brain activation in the verb generation task before 
reading skills have become automatic (as early as 7 years old) with 
reading scores at the age of 1 7 years. In this longitudinal study, chil- 
dren performed the fMRI verb generation task three times: before 
reading is fully acquired (7 years old, age range 5-8 years = T1), 
during acquisition of reading skills ( 1 1 years old, age range 9- 
14years = T2) and after reading was mastered (17 years old, age 
range 15-19 years = T3). In typical language development, read- 
ing proficiency was found to be affected by the age of reading 
acquisition, an effect that was found to decrease with age (Zevin 
and Seidenberg, 2002). Also, a greater variety in strategies used for 
reading was found in younger individuals as compared to profi- 
cient older readers (Waters et al, 1984). We therefore examined 
the relationship between fMRI activation during verb genera- 
tion, and reading scores achieved toward adulthood, at the age 
of 17 years (T3). We hypothesized that the verb generation task 
would show typical frontal, temporal, and occipital activation, 
with a trend toward left-lateralized activation with age. We also 
postulated that reading proficiency would be associated with left- 
lateralized activation in regions related to executive functions and 
reading; specifically in the frontal and occipital regions as reading 
is mastered (i.e., with development) (as was observed in Purcell 
et al., 2011). Lastly, we hypothesized that at Tl and T2 reading 
proficiency should be associated with greater activation (during 
verb generation) of regions in the cingulo- operculum network, 
because in younger children, maintaining focus on task goals 
and avoiding errors are of primary importance as reading skills 
are being acquired. In adolescence, we hypothesized that read- 
ing proficiency should be associated with greater activation of the 
fronto-parietal network at T3 after reading is mastered because 
in mature, proficient readers, rapid adaptive control mechanisms 
have been developed and support reading. 

MATERIALS AND METHODS 
PARTICIPANTS 

Sixteen children (eight boys, eight girls) were included in this lon- 
gitudinal study. Children initially entered the study at the age of 
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5-8 years and then returned annually for brain imaging and neu- 
rocognitive testing. In the present analysis, we focus on imaging 
data acquired at three time points in the 10-year longitudinal study, 
when children were of amean age of 7±0.1.03 (Tl), 1 1.53 ± 1.59 
(T2), and 17.18 ± 1.27 (T3) years. 

All the participants (N = 16) completed the verb generation 
task during fMRI at Tl and T2. In the last session (T3), only 15 
participants completed the scan due to braces in 1 of the par- 
ticipants. However, reading measures were acquired from all 16 
individuals at T3. 

The study was approved by the Institutional Review Board. 
Informed consent was obtained from the child's parent or guardian 
and assent also obtained from subjects 8 years and older. Exclu- 
sion criteria were previous neurological illness, learning disability, 
head trauma with loss of consciousness, current or past use of 
psycho-stimulant medication, pregnancy, birth at 37 weeks gesta- 
tional age or earlier, or abnormal findings at a routine neurological 
examination performed by an experienced pediatric neurologist. 
All participants were part of a parent study investigating normal 
language development in children and were considered "healthy" 
based on neurological, psychological, and structural measures 
(Holland et al, 2007; Szaflarski et al, 2012). 

All participants were native English speakers. Fifteen were right- 
handed and one was left-handed according to the Edinburgh 
Handedness Inventory (Oldfield, 1971). All participants were pre- 
screened for any conditions, which would prevent an MRI scan 
from being acquired safely. Intelligence was measured using the 
age appropriate Wechsler intelligence scales upon entry to the 
study and again 2 and 4 years later as follows: Wechsler Preschool 
and Primary Scale of Intelligence (WPPSI-R, ages below 6 years) 
(Wechsler, 1989), Wechsler Intelligence Scale for Children - Third 
Edition (WISC-III, age 6-16 years) (Wechsler, 1991), and Wechsler 
Adult Intelligence Scale - Third Edition (WAIS-III, ages 17 years 
and above) (Wechsler, 1997). Reading was assessed at the third 
session (T3 at age 17 years) using the Woodcock- Johnson Letter- 
Word reading test (Woodcock et al, 2001). In this task, the par- 
ticipants were instructed to read as accurately as possible, a list of 
words written in English, increasing in degree of difficulty. The task 
was stopped by the administrator after six reading errors in a row. 

fMRI PARADIGM 

The fMRI paradigm consisted of a covert verb generation task as 
previously detailed (Holland et al, 2001, 2007) ] using a 30-s on-off 
block design. All stimuli were presented using MacStim (White Ant 
Software, Melbourne, VIC, Australia). Stimuli were presented at a 
rate of one noun every 5 s, for six stimuli during each 30-s epoch. 
During the active epochs, the participants were asked to think of 
appropriate verbs such as "throw" or "kick" to aurally presented 
concrete nouns such as "ball." Typically children can think of two 
or three verbs associated with each noun during each response 
interval. They are instructed not to move the lips or mouth while 
covertly generating the verb responses. During the control epochs, 
participants were asked to bilaterally tap their fingers when they 
heard a modulated tone. The bilateral finger tapping task was cho- 
sen as the control for verb generation since the sensory-motor 
response to bilateral finger tapping as reported previously (Hol- 
land et al., 2001; Szaflarski et al, 2006b). In order to compare the 



control task to the auditory stimulation and response initiation 
on the verb generation task, the participants were instructed to 
sequentially tap the fingers to the thumbs on both hands simulta- 
neously when they hear a tone, in a self-paced manner. Participants 
were asked to stop tapping after touching each finger to the thumb 
twice. This control task accomplished five objectives in a pediatric 
fMRI experiment. First the auditory cue using a tone to pace the 
finger tapping controls for the auditory stimulation in the verb 
generation task. Second, tapping fingers provided activation of 
the motor strip as reference data for each participant as well as 
a method for examining differences in developmental aspects of 
BOLD activation from the neurocognitive language paradigm vs. 
a motor paradigm that should not have very strong developmen- 
tal influence over the age span tested (Schapiro et al., 2004). It 
also provided participants with a task to shift their attention away 
from generating verbs into the control epoch. In addition, paced 
sequential finger tapping requires motor planning and execution 
and serves as a control for these processes underlying verb gener- 
ation task. Finally, this control task is used to provide an indirect 
measure of participants' compliance inside the scanner. Yuan et al. 
(2009) also explored the differences in head motion children dur- 
ing fMRI while performing finger tapping vs. verb generation and 
did not find a significant difference in motion for these two phases 
of the task: a key consideration for obtaining high quality fMRI in 
pediatric subjects. The difficulty level of the fMRI task was selected 
such that children as young as 5 years old would be readily able to 
perform the task. 

In order to verify that the participants were attentive to the pre- 
sented nouns during the scan, they were given a yes/no recognition 
test involving the 25 nouns and 25 distracters. This quiz consisted 
of a sheet with a list of nouns on it: nouns that were in the task, 
and foils that were not. Participants were told to mark the words 
that they remembered hearing during the verb generation task in 
the scanner. 

Imaging 

An MRI-compatible audio visual system was used for presen- 
tation of the stimuli. Details of the techniques used to obtain 
fMRI data from younger children, as well as the success rates, 
are given in Byars et al. (2002). EPI-fMRI scan parameters 
were TR/TE = 3000/38 ms, 125 kHz, FOV = 25.6 cm x 25.6 cm, 
matrix = 64 x 64, and slice thickness = 5 mm. Twenty-four slices 
were acquired, covering the entire cerebrum. One hundred ten 
whole-brain volumes were acquired (the first 10 were discarded 
during post-processing to insure image contrast at relaxation equi- 
librium) for a total scan time of 5 min 30 s. Techniques detailed 
elsewhere (Byars et al., 2002) were used to acclimatize the par- 
ticipants to the MRI procedure and render them comfortable 
inside the scanner. Soft head restraints were used to minimize head 
motion. In addition to the fMRI scans, whole-brain Tl weighted 
MP-RAGE scans were acquired for anatomical co-registration. 
All imaging was performed using a 3 T, head only MRI scanner 
(Bruker Medspec 30/60). 

DATA ANALYSIS 

Data were analyzed using Cincinnati Children's Hospital Image 
Processing Software (CCHIPS), written in IDL (Research Systems 



Frontiers in Human Neuroscience 



www.frontiersin.org 



June 2014 | Volume 8 | Article 447 | 3 



Horowitz-Kraus et al. 



Neural-correlates to verb-generation and reading 



Inc., Boulder, CO, USA). Image data were corrected for Nyquist 
ghosts and geometric distortion using multi-echo reference 
method (Schmithorst et al., 2001), and was motion-corrected 
using pyramid co-registration (Thevenaz et al., 1998): we per- 
formed three-dimensional ridged body transformation to align 
the volumes. This resulted in six motion parameters. These para- 
meters were included as regressors in the first-level general linear 
model (GLM) analysis for each subject's data. To ensure that the 
level of motion did not differ across the three sessions, the mean 
motion was reported as well as included in an analysis of vari- 
ance (ANOVA). We also used a within group f-test analysis to 
compare the motion levels between the task (when generating 
verbs) and the contrast conditions (when tapping fingers). In 
addition, time points with excessive motion were rejected from 
the post-processing pipeline. We used a mutual information cost 
function for rejecting motion corrupted frames of fMRI data as 
previously described for this same cohort of longitudinal subjects 
(Szaflarski et al, 2006). As mentioned in our previous studies, 
covert and overt verb generation resulted in the same head motion 
in this same cohort of children (Vannest et al., 2009). Although the 
data included in the current analysis are from a slightly different 
block-periodic version of the verb generation task, we did not find 
differences in head motion during the covert verb generation and 
finger tapping conditions (Yuan et al., 2009). All data met the cri- 
terion of median voxel displacement in the center of the brain 
<2 mm. The fMRI data were transformed into stereotaxic space 
(Talairach and Tournoux, 1988) using a linear affine transforma- 
tion (see Muzik et al., 2000). The use of the Talairach standard for 
children aged 5 years and above has been shown to produce mini- 
mal errors in co-registration for group analysis. Consequently, we 
are confident that using the adult Talairach standard in the present 
study will not substantially alter the findings of the current study 
(Wilke et al, 2002; Altaye et al, 2008). Note that the imaging data 
for this study were acquired over a period of 12 years and that the 
methods for fMRI analysis have improved considerably over this 
period of time. After data acquisition on all participants was com- 
pleted, 12 years after the longitudinal study began we processed all 
of the data from all participants specifically for this report using 
the analysis pipeline as described above. 

GROUP ACTIVATION MAP 

A GLM and random-effects analysis were used in order to detect 
the group activation for Tl, T2, and T3 children for the contrast 
between conditions of verb generation > finger tapping contrast. 
Images of the f-maps generated by this first-level contrast were 
thresholded for visual inspection purposes to p < 0.001, corrected 
via Monte Carlo simulation using a cluster of 30 voxels (Forman 
et al, 1995). 

Regression 

We used regression in second level analyses in order to examine 
the association between reading scores from the "Letter-Word" 
subtest (from the Woodcock-Johnson battery) from the third 
session (T3) and fMRI activation during the verb generation 
task at T1-T3. Following the first-level computation of individ- 
ual (--score maps for the verb generation > finger tapping con- 
trast, whole-brain voxel-wise group maps, were computed using 



Spearman's rank-correlation coefficient. These maps were based 
on unthresholded f-maps from the first-level analysis and each 
subject's behavioral score included as a predictor of activation (e.g., 
"Letter- Word" subtest from the Woodcock-Johnson III). Images 
were thresholded to p < 0.001, uncorrected. Cluster size for each 
group was 35 voxels. 

Lateralization index calculation 

A lateralization index (LI) was also calculated for each subject in 
the three age groups based on the unthresholded individual z- 
score maps. The LI represents a relative hemispheric difference for 
an individual that is self-normalizing in terms of relative BOLD 
activity. Since this index is non-dimensional, it provides a con- 
venient means of comparing groups of subjects across different 
scanning sessions. As noted later, the brain areas included for LI 
calculation were determined so that they were consistent across 
the groups while still conforming to the traditional standard of 
regions of interest (ROI) selection for language lateralization. 
Every effort was made to insure that the methodologies applied 
in the LI calculation of the three subject groups were compatible. 

Several ROIs were chosen for the purpose of the current analysis 
using an LI toolbox implemented in CCHIPS. Language-network 
regions were defined for LI calculation based on the composite 
map of activation obtained from the fMRI verb generation task 
from all subjects at T3. First, we chose a frontal ROI correspond- 
ing to the traditional "Broca's area"; we retained all active voxels 
within and contiguous to this area, so the frontal ROI included the 
entire IFG (BA 44, 45, 47) and extensions into additional dorso- 
lateral prefrontal regions, including portions of the MFG (BA 46, 
48, 49) and precentral gyrus (BA 46). We also selected a temporal 
language ROI based on the traditional "Wernicke's area," defined 
functionally by the group activation map. This ROI extended from 
the temporal plane (BA 41, 42) inferiorly through the middle tem- 
poral gyrus (MTG) to the margin of the ITG (including portions 
of BA 22, 21, 37) but did not include posterior aspects (BA 22, 
39). ROIs for the right hemisphere homologs were established by 
reflecting the coordinates to the right hemisphere. 

Using the same approach, we defined ROIs for reading 
(fusiform gyrus, BA 37) and for executive functions (dorsolat- 
eral PFC, BA 9; and the anterior PFC, BA 10) (after Purcell et al, 
2011). We retained all active voxels within and contiguous to these 
areas and defined them as two additional separate ROIs on the left 
and for the right hemisphere homologs based on the correlation 
map of activation obtained from the fMRI verb generation task 
correlated with the WJ-III reading scores from all subjects at T3. 

Lateralization indexes for language-network regions were com- 
puted both for the composite f-maps and for the regression maps; 
the Lis for the reading and executive functions-related regions were 
then computed separately from the regression maps. Only voxels 
with z-scores greater than or equal to the mean z-score within an 
ROI for each individual subject were used in the calculation of Lis. 
Voxels above this mean z-score threshold were counted and an LI 
was defined as the difference in the number of activated voxels, 
summed independently for the left and right ROIs, divided by the 
sum total of active voxels in the left plus right ROI. According 
to this formula, a positive LI indicated left hemisphere lateraliza- 
tion and a negative number indicated right lateralization. These 
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LI values were used to construct the lateralization growth curves 
(a linear regression line). 

RESULTS 

NEUROPSYCHOLOGICAL TESTING 

Neuropsychological testing results verified that the par- 
ticipants had normal to above average scores in read- 
ing and verbal abilities. Mean Wechsler verbal and non- 
verbal IQ remained stable during development. At Tl: ver- 
bal scaled IQ= 115.9± 9.23 (range = 104-139), non-verbal 
scaled IQ=112.9± 12.23 (range = 96-130), at T3: verbal 
scaled IQ= 109.75 ±0.37 (range = 99-139), non-verbal scaled 
IQ= 114.91 ±14.71 (range = 95-123). 

VERB GENERATION - POST-SCANNING MEASUREMENT 

The average number of correct responses for the verb genera- 
tion post-test outside the scanner at T2 and T3 were 12.4/25 ± 3.6 
and 15.05/25 ± 2.4, respectively. Range of possible scores is 0-25 
correct choices. 

RANDOM-EFFECTS ANALYSIS 

The statistical parameter maps for verb generation > finger tap- 
ping for T1-T3 in the longitudinal cohort shown in Figures 1-3 
are consistent with previous studies using this task (Holland et al, 
2007; Szaflarski et al., 2012). While there are many areas of acti- 
vation that are common across ages in the cohort, there are also 
variations. 

77 children 

Significant activation was found in the left and right insula (BA 
13) and left IFG (BA 13, 9), left fusiform gyrus (BA 20), left STG 
(BA 22), left MTG (BA 21), left superior frontal gyrus (SFG) 
(BA 6), left cingulate gyrus (BA 32), left thalamus, right pre- 
central gyrus (BA 6), and left and right parahippocampal gyrus. 
Talairach coordinates of cluster centroids are listed in Table 1 (see 
Figure 1). 

72 children 

Significant activation was observed in the left and right MFG (BA 
46, 9, 10), left IFG (BA 47), left MTG (BA 21, 22), right IFG (BA 




FIGURE 1 | Composite fMRI activation maps for the verb generation 
(verb generation > finger tapping) in Tl children ( N = 1 6). The contrast is 
significant at p < 0.001 (corrected), slice thickness is 5 mm for these 
contiguous slices. Slices range from z= 13 to 24 in the Talairach frame. 
Cluster size is 30 voxels. Higher significance is indicated in hotter color (t 
threshold ranged from 3 to 6, see scale in the bottom of the figure). Figure 
is presented in a radiological orientation (L= R, R = L). 



13), and left anterior cingulate (BA 32). Talairach coordinates of 
cluster centroids are listed in Table 1 (see Figure 2). 

73 children 

Significant activation was observed in the left STG (BA 22), left 
MFG (BA 46, 9, 8), left parahippocampal gyrus (BA 36), left (BA 
13) and right IFG (BA 47), and right cuneus (BA 30). Talairach 
coordinates of cluster centroids are listed in Table 1 (see Figure 3). 

REGRESSION ANALYSIS MAPS 

Regression analyses maps using the "Letter- Word" standard scores 
from the Woodcock-Iohnson III at T3 as a predictor, revealed 
significant positive correlations with the contrast of verb gen- 
eration > finger tapping in the children at three different age 
points. 

T1 children 

A significant positive correlation was observed in the left insula 
(BA 13), left and right MFG (BA 8, 6), left posterior cingulate (BA 
30), left anterior cingulate (BA 24), left thalamus, right precentral 
gyrus (BA 44), right superior temporal lobe (BA 41), right MTG 
(BA 37), right cuneus (BA 7), and right posterior cingulate (BA 23, 
30). Talairach coordinates of cluster centroids are listed in Table 2 
(see Figure 4). 




FIGURE 2 | Composite fMRI activation maps for the verb generation 
(verb generation > finger tapping) in T2 children {N = 1 6). The contrast is 
significant at p < 0.001 (corrected), slice thickness is 5 mm for these 
contiguous slices. Slices range from z= 13 to 24 in theTalairach frame. 
Cluster size is 30 voxels. Higher significance is indicated in hotter color (f 
threshold ranged from 3 to 6). Figure is presented in a radiologica 
orientation (L= R, R = L). 




FIGURE 3 | Composite fMRI activation maps for the verb generation 
(verb generation > finger tapping) in T3 children {N = 1 5). The contrast is 
significant at p < 0.001 (corrected), slice thickness is 5 mm for these 
contiguous slices. Slices range from z= 13 to 24 in theTalairach frame. 
Cluster size is 30 voxels. Higher significance is indicated in hotter color (f 
threshold ranged from 3 to 6). Figure is presented in a radiologica 
orientation (L = R, R = L). 
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Table 1 | Group composite contrast (verb generation > finger tapping) inT1,T2, andT3 participants. 



Region BA Cluster size Cluster centroid 



X Y Z 
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I oft mfonor frontal n\/n ic 
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13 


66 
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27 


14 




9 


150 


-39 


8 


29 


Left superior temporal gurus 


22 


155 


-51 


-24 


3 


Left middle temporal gyrus (sub-cluster) 


21 


26 


-48 


-36 


-5 


Left fusiform gyrus (sub-cluster) 


20 


22 


-44 


-37 


-10 


Left superior frontal gyrus 


6 


32 


-8 


15 


57 


Left parahippocampal gyrus 




202 


-26 


2 


-15 


left thalamus 




142 


-8 


15 


57 


Left cingulate gyrus 


32 


64 


-7 


22 


39 


Right precentral gyrus (frontal lobe) 


6 


41 


42 


-7 


30 


Right parahippocampal gyrus 




224 


26 


4 


-18 


Right insula 


13 


48 


42 


-18 


13 




Left middle frontal gyrus 


46 


15 


-38 


46 


15 


Left middle frontal gyrus (sub-cluster) 


9 


63 


-43 


12 


27 


Left inferior frontal gyrus (sub-cluster) 


47 


88 


-40 


18 


-13 




13 


345 


-44 


23 


7 


Left middle temporal gyrus 


22 


94 


-60 


-35 


5 


Left middle temporal gyrus (sub-cluster) 


21 


21 


-58 


-35 


-13 


Left anterior cingulate 


32 


91 


-2 


19 


40 


Right inferior frontal gyrus 


13 


69 


38 


13 


-13 


Right middle frontal gyrus 


10 


52 


40 


38 


22 


Left parahippocampal gyrus (limbic lobe) 


36 


82 


-39 


-37 


-6 


Left superior temporal gyrus 


22 


87 


-47 


-22 


4 


Left inferior frontal gyrus 


13 


284 


-38 


24 


8 


Left middle frontal gyrus (sub-cluster) 


9 


78 


-38 


14 


26 




6 


28 


-37 


4 


40 


Right middle frontal gyrus 


8 


158 


0 


19 


43 


Right cuneus 


30 


56 


1 


-68 


7 


Right Inferior frontal gyrus 


47 


44 


37 


23 


1 



12 children 

A significant positive correlation was observed in left and right 
fusiform gyri (BA 37, 19), left STG (BA 41), left anterior cingu- 
late (BA 32), and the MFG (BA 10). In the right side, a significant 
correlation was observed in the caudate, hypothalamus, and hip- 
pocampus. Talairach coordinates of cluster centroids are listed in 
Table 2 (see Figure 5). 

73 children 

A significant positive correlation was observed in the left STG (BA 
41, 38), left MFG (BA 9, 10), left insula (BA 13), left fusiform 
gyrus (BA 20), left parahippocampal gyrus (BA 34), and the 
left precuneus (BA 7) and cuneus (BA 18). A significant cor- 
relation was also observed in the right MFG (BA 6), precuneus 
(BA 31), posterior cingulate (BA 30), and lingual gyrus (BA 18). 
Talairach coordinates of cluster centroids are listed in Table 2 (see 
Figure 6). 



AGE-RELATED CHANGES IN LATERALIZATION 

In Figure 7 the LI for each group is plotted as a function of age 
(for Tl, T2, and T3) for frontal and temporal language-network 
ROIs illustrated in the composite maps in Figures 1-3. A signif- 
icant age-related increased left-lateralization was found for the 
frontal ROI (i.e., "Broca") (R 2 = 0.92). Although we observed a 
leftward lateralization trend in the temporal region as well (i.e., 
"Wernicke"), age-related change in this posterior ROI was more 
modest (_R 2 = 0.02) . Results suggest that age accounted for as much 
as 92% of individual variance in LI in the activation in Broca and 
for 2% of individual variance N LI in the activation in Wernicke. 
See upper part of Figure 7. 

An examination of the language ROIs from the regression 
maps with reading scores in Figures 4-6, revealed increased left- 
lateralization along development as shown in the lower panel of 
Figure 7 for the frontal (green triangles, _R 2 = 0.99) and tem- 
poral (red squares, R 2 = 0.9) language-network ROIs. Greater 
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Table 2 | Regression analysis of MRI data from the verb generation task atT1,T2, andT3 participants with Letter-Word reading standard scores 
atT3 



Region BA Cluster size Cluster centroid 



X Y Z 
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Right cuneus (occipital lobe) 


7 


51 


9 


-66 


30 


Left fusiform gyrus (temporal lobe) 


37 


48 


-43 


-50 


-10 


Left superior temporal gyrus (sub-cluster) 


41 


20 


-38 


-42 


5 


Left anterior cingulate 


32 


47 


-5 


31 


20 


Left middle frontal gyrus 


10 


36 


-9 


50 


5 


Left fusiform gyrus (occipital lobe) 


19 


40 


-38 


-73 


-12 


Right caudate (temporal lobe) 




25 


36 


-30 


-5 


Right culmen (anterior) (sub-cluster) 




25 


25 


-29 


-25 


Right hypothalamus 




37 


1 


-2 


-10 


Right hippocampus (temporal lobe) (sub-cluster) 




17 


34 


-29 


-10 


Right fusiform gyrus (occipital lobe) 


19 


35 


34 


-77 


-12 




Left superior temporal gyrus (temporal lobe) 


38 


44 


-39 


14 


-22 


Left fusiform gyrus (temporal Lobe) 


20 


96 


-36 


-38 


-15 


Left parahippocampal gyrus (limbic lobe) 


34 


41 


-16 


-7 


-18 


Left superior temporal gyrus 


41 


111 


-44 


-37 


4 


Left insula 


13 


100 


-37 


-20 


9 


Left insula (sub-cluster) 


13 


44 


-35 


-13 


14 


Left middle frontal gyrus (frontal lobe) 


10 


86 


-37 


37 


9 


Left middle frontal gyrus 


9 


97 


-33 


15 


27 


Left cuneus (occipital lobe) 


18 


58 


-6 


-76 


23 


Left precuneus (parietal lobe, sub-cluster) 


7 


16 


-19 


-55 


56 


Right medial frontal gyrus (frontal lobe) 


6 


81 


6 


-19 


53 


Right precuneus 


31 


76 


18 


-74 


28 


Right posterior cingulate 


30 


75 


31 


-74 


12 


Right lingual gyrus (occipital lobe) 


18 


72 


17 


-54 


5 



age-related changes in lateralization, were found for the reading 
regions as well (turquoise X, fusiform gyrus; _R 2 = 0.97), and for 
the executive functions regions (purple X, BA 9: R 2 = 0.83, BA 
10: R 2 = 0.96). See lower part of Figure 7. Results suggest that 
age and reading proficiency accounted for 99% of individual vari- 
ance in LI in the activation in Broca and for 90% of individual 
variance in LI in the activation in Wernicke. It also emphasized 



that age and reading proficiency accounted for 97% of individual 
variance in LI in the activation in the Fusiform gyrus and for 83 
and 96% of individual variance in LI in the activation in executive 
functions regions (BA 9 and 10, respectively). These findings indi- 
cate that those regions positively correlated with reading scores 
became more left-lateralized with age. This was true for function- 
ally defined regions of significant correlation for reading scores 
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FIGURE 4 | Regression maps showing a positive correlation for 
activation during verb generation scans inT1 with the Letter-Word 
score from the Woodcock-Johnson III [N= 16) inT3. All activated pixels 
meet significance threshold of p < 0.05, corrected. Slices range from z = 7 
to 24 in theTalairach frame. Cluster size is 35 voxels. Higher significance is 
indicated in hotter color (r value ranged from 0.2 to 0.5, see scale in the 
bottom of the figure). Figure is presented in a radiological orientation (L= R, 
R = L). 




FIGURE 5 | Regression maps showing a positive correlation for 
activation during verb generation scans atT2 and the Letter-Word 
score from the Woodcock-Johnson III [N= 16) fromT3. All activated 
pixels meet significance threshold of p < 001 uncorrected. Slices range 
from 2 = 1 to 24 in theTalairach frame. Cluster size is 35 voxels. Higher 
significance is indicated in hotter color (r value ranged from 0.2 to 0.5). 
Figure is presented in a radiological orientation (L = R, R = L). 



and brain activation as well as for regions that were active in the 
first level analysis in Broca and Wernicke's areas. The fact that this 
lateralization trend is similar for all of these regions (slope of LI vs. 
time in the range of 0.37-0.59) suggests that these regions are all 
connected in support of reading and language skill and are devel- 
oping at a similar pace, even though the regression coefficient is 
not significant for the primary regions. 

DISCUSSION 

In the current study, we aim to examine the relationship between 
reading proficiency and the maturation of neural networks related 
to executive functions and language, using the verb generation 
task. The verb generation task draws on language and executive 
skills that improve during development (Missier and Crescentini, 
2011) and correlate with reading ability (Snowling et al, 1997); 
thus, we were able to study participants in the youngest group 




FIGURE 6 | Regression maps showing a positive correlation for 
activation during verb generation scans fromT3 and the Letter-Word 
score from the Woodcock-Johnson III [N= 15) atT3.The activated pixels 
meet significance at p < 0.001 , uncorrected, higher significance is indicated 
by hotter color. Slices range from z = 7 to 24 in theTalairach frame. Cluster 
size is 35 voxels. Higher significance is indicated in hotter color (r value 
ranged from 0.2 to 0.5). Figure is presented in a radiological orientation 
(L=R, R=L). 



(Tl) who entered the study before they formally acquired reading 
skills or before reading was fully acquired, and examine activa- 
tion in regions supporting language and executive skills without 
requiring a reading task. 

In line with our hypothesis, our results demonstrate a trend 
of more left-lateralized activation during the verb generation task 
from Tl to T3. In other words, age was found to account for 
the increasing left LI in language-related regions. Specifically, we 
observed increasingly leftward lateralization of the inferior/middle 
frontal region with development, as reported previously in a super- 
set including the participants of the current analysis (Holland 
et al., 2001, 2007). This trend was less robust in the temporal 
language region, as previously observed, due to a weaker reliance 
on auditory-phonological-related regions during this task (Hol- 
land et al., 2007). A similar trend of left-lateralized activation 
with development was observed by Gaillard et al. (2000) while 
administering the verb generation task to participants age 10.7 
and 28.7 years. This developmental shift to the left hemisphere 
for linguistic tasks was also reported in other studies (Szaflarski 
et al., 2006a) and reflects the specialization of the left hemisphere 
for language. These results are not confounded by differences in 
motion across age or between the different task conditions (i.e., 
generating verbs vs. tapping fingers). Although the head motion 
data were regressed out of the analysis, the absence of the effect 
of age on motion level, confirms that these findings of increased 
left-lateralization with age are not because of motion differences 
(see also Gaillard et al, 2003). 

The main goal of the current study was to identify the develop- 
mental changes in the activation of the neural circuitry supporting 
language and executive functions (during the verb generation 
task in Tl, T2, and T3) corresponding to proficient reading at 
T3. Our results confirm our second hypothesis: reading profi- 
ciency was found to be associated with left-lateralized activation 
in regions related to language and profound left-lateralization 
in executive functions (BA 9, 10) and reading (BA 37). In gen- 
eral, greater reading scores in Tl were associated with a bilateral 
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Lateralization index for Broca and Wernicke regions during the verb generation task 

in T1, T2andT3 




Lateralization index for language (Broca, Wernicke), frontal (BA 9 and 10) and occipital (BA 37) regions 
in a verb generation fMRI data correlated with reading scores at T3 




FIGURE 7 | Lateralization index (LI) for the composite (upper) 
and the regression maps (lower). LI was calculated on 
language-related regions (Broca and Wernicke) in the composite 
maps in T1— T3. LI was also calculated on language (Broca and 



and more widespread brain activity, in addition to the "classical" 
reading regions. This involved the right STG and MTG, which 
are related to phonology and bilateral activation of production 
and semantic-related regions such as the insula. These findings 
extend our previous results (Holland et al., 2001, 2007) by show- 
ing that age and reading proficiency account for increasing left- 
lateralization in regions related to language, reading, and executive 
functions. This is in line with findings suggesting lesser reliance 
on phonological and semantic processes and greater reliance on 
orthographic processes during reading proficiency (Ehri, 2005), 
and neuroimaging findings of greater left-lateralization during 
reading along development (Turkeltaub et al., 2003). With age, 
the association of greater reading scores and brain activation in 
the verb generation task was found to be more specific to execu- 
tive functions and reading-related regions. At T2 and T3, greater 



Wernicke), executive functions (dorsolateral prefrontal cortex, 
anterior prefrontal cortex), and reading regions (fusiform gyrus). 
Standard deviations representing the distribution of individual Lis 
are indicated. 



reading scores were associated with more leftward LI for language- 
related regions (frontal and temporal) than at Tl, but even more 
leftward LI scores for frontal and occipital regions related to exec- 
utive functions and reading, respectively. Since better reading is 
associated with greater reliance on word recognition (i.e., "ortho- 
graphic route") and less on decoding (i.e., "phonological route"), 
these results might demonstrate the decreased role of language 
regions (i.e., phonology-auditory route) and the increased role of 
the orthography-visual route in reading development (for more 
information see Coltheart et al., 2001; dual route model). 

Positive correlation of word reading scores and the activation 
of the fusiform gyrus was observed only in T2 and T3, whereas 
children in the Tl group showed greater right cuneus and left 
precuneus activation. We suggest that because younger age chil- 
dren have not reached reading proficiency and the mental lexicon 
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Table 3 | Motion measures, mean, and standard deviation (SD) for 
children atT1,T2, andT3. 

Mean (SD) Mean (SD) Mean (SD) (Test 

(total: generating (generating (tapping (generating 

verbs and verbs) fingers) verbs vs. 

tapping fingers) tapping 

fingers) (ns) 

T1 children 0.3(0.22) 0.33(0.22) 0.3(0.21) 0.38 

T2 children 0.39 (0.36) 0.43 (0.53) 0.37 (0.22) 0.26 

T3 children 0.33(0.3) 0.32(0.31) 0.34(0.3) 0.43 



Means and standard deviations (SDs) for motion measures for the total task (gen- 
erating verbs and finger tapping!, the task condition (generating verbs), and for the 
contrast condition (finger tapping! across the different age groups. No significant 
differences were found between the different age groups for the total motion 
parameter as well as between the task and contrast conditions (left column). 

has not fully been established, there is heavier load on regions 
related to phonological processing (e.g., STG, MTG) and less on 
regions related specifically to automatic word recognition (e.g., 
fusiform gyrus) earlier in development. At T2, children are using 
orthographic processing regions as they become more proficient 
readers (Ehri, 2005) and therefore show a correlation between 
reading proficiency with bilateral fusiform gyri activation. Finally, 
left sided fusiform gyrus activation is observed at T3, as expected 
from proficient and automatic readers (van der Mark et al., 201 1). 
The activation of the fusiform gyrus, lingual gyrus, and other 
occipital regions may be related to visual imagery during the 
verb generation task (see also Karunanayaka et al, 2011). This 
process was also found to correspond to greater reading proficiency 
(Horowitz-Kraus et al, 2013). 

Note that the ANOVA we performed to assess the effect 
of age on motion did not reveal significant age differences 
[F (3,45) = 0.48, p > 0.05]. Namely, the level of motion did not 
differ across the age points. No significant differences were found 
within groups either, suggesting that motion level for the task and 
contrast conditions was relatively equal (see Table 3). 

The frontal lobe matures with age, along with concomitant 
improvements in executive functions (Segalowitz and Davies, 
2004). Giedd observed that gray and white matter maturation, 
which includes completion of myelination and synaptic pruning, 
in the parietal and frontal regions peaks at 16 years old, so the 
cognitive abilities centered in these regions are mature as well 
(Giedd et al, 1999). Some later findings point at an inverse rela- 
tionship between cortical gray matter density reduction and brain 
growth primarily in the superior frontal regions that control exec- 
utive functioning even in 30-year-old adults (Sowell et al., 2001). 
Executive functions, which support inhibition, WM, planning, and 
attention, likely develop throughout adolescence together with the 
maturation of the frontal and parietal cortices (Casey et al, 2000; 
Luna etal.,2010). 

Since our sample is children in the age range of 7-17 years, 
frontal lobe maturation is still in progress, and verb genera- 
tion may demand a reliance on different neural-circuits related 
to executive functions for Tl and T2 children as compared to 



T3. We hypothesized that the fronto-parietal network may be 
more engaged in older children, and our finding of greater left 
frontal and parietal (precuneus) activation correlated with read- 
ing proficiency at T3 supports this hypothesis. As suggested by 
the dual-networks model, the fronto-parietal network may be in 
charge of rapid adaptive control and attention allocation, which 
are part of executive functions (Dosenbach et al., 2008). This 
network has been previously found to correlate with higher read- 
ing scores (Vogel et al., 2014). As adolescents become proficient 
readers, rapid monitoring is required to process the written stim- 
uli in order to achieve fluent reading. This same rapid control 
mechanism that contributes to fluent reading is also engaged to a 
greater degree during verb generation in those participants who 
are the most proficient readers. The cingulo-opercular network, 
supporting maintenance of task goals and error monitoring, is also 
observed, in part, at T3: the insula was also positively correlated 
with higher reading scores. 

We also found that the cingulo-opercular network was more 
engaged in younger children; there was activation of the insula 
at Tl and of the anterior cingulate cortex at Tl and T2. These 
results suggest that greater task maintenance and error monitoring 
is needed by younger children - as they perform the verb gener- 
ation task, semantic associations among words are not as well- 
established, and selecting appropriate responses involves greater 
resources. This is consistent with a previous finding from the verb 
generation task when nouns presented to the participants involved 
greater conflict and the need for more attentional resources to 
resolve this conflict (Barch et al, 2000). 

The activation of the anterior cingulate cortex in Tl and T2 
was also positively correlated with higher reading scores at the 
age of 17 years (T3). This suggests that children who go on to 
become more proficient readers engage this region to a greater 
degree even early in development. Before reading becomes auto- 
matic, there are several competing representations of new words 
as they are encountered, until the orthographic route becomes sta- 
ble (Coltheart et al, 2001; dual route model). Young children who 
devote increased resources to resolving these conflicting repre- 
sentations may ultimately become better readers, suggesting that, 
more generally, executive functions play a role in the reading acqui- 
sition process. We suggest that by adolescence, as linguistic and 
executive abilities increase, conflict may decrease, resulting in the 
absence of activation in the anterior cingulate cortex in the verb 
generation task at T3. The anterior cingulate cortex also sup- 
ports error detection and is active when making reading errors 
(Horowitz-Kraus and Breznitz, 2008), and the absence of associ- 
ation between reading ability and the activation in this region at 
T3 years might be due to less error detection activation during the 
verb generation task as well as in reading and less need for set- 
maintenance during these processes, which become automatic. 
Vogel et al. (2013) found no developmental changes in fronto- 
parietal and cingulo-opercular networks when examining these 
networks using resting-state functional connectivity. Our regres- 
sion maps as well as the LI values from the frontal ROIs suggest 
that while these executive functions networks remain equally active 
across development, when examining a specific task; one net- 
work may be more active than another at different developmental 
stages. 
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To conclude, our results support the utilization of language 
and executive functions regions for development of proficient 
reading. These domains also have a differential role of execu- 
tive functions in reading development: greater role of the set- 
maintenance network before reading becomes automatic (Tl and 
T2) and greater activation of the rapid, adaptive control net- 
work when reading becomes proficient. Our results strengthen 
previous studies, which pointed at the crucial role of executive 
functions in the reading process, highlighting the differential need 
in sub-components of executive functions along development. 
Clinically, the results provide neuro-functional support for the 
importance of intact executive functions in future reading devel- 
opment, by means of early diagnosis and intervention, especially 
among populations with reading impairments (dyslexia, atten- 
tion deficit hyperactive disorder). These conclusions should be 
taken in the context of the current study's limitations: (1) We 
did not acquire fMRI data during a reading task, so we can- 
not compare the regression maps with reading scores and verb 
generation to actual activation during reading in the same par- 
ticipants; (2) there is a reading exposure at 7 years old (and 
especially since the age range of this group was relatively wide) 
so the youngest group of participants is not completely "naive" 
to reading. Adding to this ambiguity in the sample is the lack of 
assessment of reading ability at the younger age points. An addi- 
tional study enrolling children younger than 7 years old should 
be done to verify whether the same phenomenon exists also 
in a younger age. (3) We did not assess the verb generation 
ability after the scan at 7 years old since the children were too 
young. (4) The only reading measure that was used in the cur- 
rent study was a word reading task. Although word reading is 
highly related to reading comprehension, these two abilities rely 
on different basic cognitive abilities (see Christopher et al., 2012 
for more information) as well as different white matter tracts 
(Horowitz-Kraus et al, 2013). Also, whereas an untimed word 
reading task, as was used in the current study involves phono- 
logical and orthographical processing, a timed word recognition 
task relies additionally on speed of processing (see Wolf and Bow- 
ers, 1999; Hart et al, 2010). Therefore, a correlation of the fMRI 
data from the verb generation task with a timed word reading 
task may result in less activation of phonological-related regions 
(e.g., STG). A future study should look at the longitudinal dif- 
ference in the fMRI verb generation correlates with these two 
measures. Despite these limitations, it is important to note that 
there are only a few studies assessing the relationship between 
executive functioning and reading and the current study is the 
first to do so longitudinally. Therefore, our findings (regional cor- 
relation of later reading abilities with early fMRI results) offer 
the possibility of a novel predictive biomarker for future reading 
ability. 
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